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Abstract. One of the aims of systems biology is to build multiple lay¬ 
ered and multiple scale models of living systems which can efficiently 
describe phenomena occurring at various level of resolution. Such mod¬ 
els should consist of layers of various microsystems interconnected by a 
network of pathways, to form a macrosystem in a consistent way; that 
is, the observable characteristics of the macrosystem should be, at least 
asymptotically, derivable by aggregation of the appropriate features of 
the microsystems forming it, and from the properties of the network. In 
this paper we consider a general macromodel describing a population 
consisting of several interacting with each other subgroups, with the 
rules of interactions given by a system of ordinary differential equations, 
and we construct two different micromodels whose aggregated dynamics 
is approximately the same as that of the original macromodel. The mi¬ 
cromodels offer a more detailed description of the original macromodel’s 
dynamics by considering an internal structure of each subgroup. Here, 
each subgroup is represented by an edge of a graph with diffusion or 
transport occurring along it, while the interactions between the edges 
are described by interface conditions at the nodes joining them. We 
prove that with an appropriate scaling of such models, roughly speak¬ 
ing, with fast diffusion or transport combined with slow exchange at the 
nodes, the solutions of the micromodels are close to the solution to the 
macromodel. 
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1. Introduction 

In the recent paper m, the authors have proposed an interesting and exciting 
paradigm of developing theoretical biology through a proper mathematization 
of systems biology. It involves an interplay of many disciplines belonging to 
mathematics, biology and to their overlap, such as the theory of evolution, 
immune competition, mutation and selection, kinetic theory, evolutionary 
game theory, multi-scale methods and networks. In the words of D. Noble, 
Ref. [ 32 ], PP- and 129, 

We are looking towards a mature theory of biological systems-level 
interactions [...]. The task of systems biology is first to unravel 
these interactions and then to develop theories to account for them, 
and so to lay bare their logical underpinnings. 

One of the important goals of integrative systems biology is to 
identify the levels at which the various functions exist and operate. 

In other words, the aim of systems biology is to build universal, hierarchical 
models of biological phenomena which would include all levels of organization 
of living matter. For instance, for malaria, one would like to have a model 
including the dynamics of the plasmodium, through the cell, tissue, individual 
and ending at the population, or even metapopulation, levels. 

1.1. Multi-scale models in systems biology. 

Systems biology has reached out to many areas of mathematics, even includ¬ 
ing disciplines such as category theory [21], or logic and set theory [33] which, 
for many years, have been regarded as belonging to pure mathematics. How¬ 
ever, modelling solely based on these disciplines results in a static picture of 
the system, which is not satisfactory for many applications. To provide the 
required dynamical features we enhance the model by incorporating differ¬ 
ential equations into the description. These equations are intended to model 
the evolution of particular building blocks of the system and typically they 
are interlinked by a complex network spanning several time and size scales. 
Usually the existence of different scales, or levels, in the model is revealed by 
the presence of small (or large) scaling parameters representing the ratios of 
the times typical for processes occurring in particular building blocks, or the 
ratios of the sizes of objects involved in them. In full generality, such models 
are too complex to allow for any robust analysis and thus it is of interest 
to be able to focus on the scale relevant to a particular aspect of the sys¬ 
tem’s behaviour by selecting an appropriate simpler sub-model. However, in 
many applied sciences often there is an expectation that complex multi-scale 
systems can be described by plug-and-play type models. In such models, fo¬ 
cusing on a required level amounts just to switching off the unwanted scales 
by setting relevant scaling parameters to 0. Unfortunately, in most cases it is 
impossible. Even in the text quoted above we see that there is a clear under¬ 
standing that all levels of organization are interconnected. In other words, it 
is recognized that there is an interdependence between various scales in the 
model so that, whenever we focus on a sub-model acting at some level, there 


Asymptotic state lumping in transport and diffusion problems on networks! 


Ex 

E3 

e 5 

E 8 

o 

k a 

o 

v 2 

0 e 7 

Ex - 

— ^3 

0 

^5 



E w \ Ee 




JC % 



Figure 1 . An example of a directed graph. 


should be a ‘shadow’ of all scales that were switched off. The problem is ex¬ 
acerbated by the fact that dynamics at different levels usually is described by 
equations which are not fully compatible with each other, especially close to 
the boundaries between the levels. In other words, setting a particular scaling 
parameter equal to zero may dramatically change the type of the equation 
and render the problem ill-posed. Thus, moving between the scales, if possi¬ 
ble at all, cannot be accomplished by simply setting the appropriate scaling 
parameter to 0 but requires complicated limit procedures which lead to the 
so-called singular limits of the involved equations. 

1.2. Multi-scale models on networks. 

The paper Ref. [TO] offers a survey of a wide range of mathematical methods 
which allow for dealing with such complex and interrelated models. In the pre¬ 
sented paper we shall focus on particular aspects of the proposed approach, 
namely on multi-scale dynamical systems related by a network and their sin¬ 
gular limits. Such systems provide a detailed description of a model at what 
we refer to as the micro-scale, while the singular limits give an aggregated 
description of the model at the macro-scale. This procedure is often referred 
to as the asymptotic state lumping. In this paper we focus on two types of 
dynamical systems, transport and diffusion and, to keep technicalities to the 
minimum, we restrict ourselves to linear problems. 

To set the stage, we begin with a brief description of how complex, 
more detailed, systems are derived from more crude building blocks. When 
we model a complex, interlinked system of subpopulations, we often begin 
by considering a static graph (or often a directed graph), where the edges 
represent connections between the subpopulations and the intensity of inter¬ 
actions are given by the weights associated with the edges, see e.g. Ref. m- 
As we mentioned above, static models usually fall short of what is expected 
from them. Therefore the next step is to assume that a subpopulation, local¬ 
ized at a particular node, changes due to interactions with subpopulations 
at the nodes connected with that node. It results in a system of ordinary 
differential equations, where, in the linear case, the weights are the rates at 
which a particular node influences the changes at the connected nodes. 
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Figure 2. A circuit realization of the Aristazabal 
Glavinovic model. 

1.2.1. First examples. 

Example 1.1. A basic mutation model. Consider the population described by 
v = (vi, ..., v m ), where Vj, 1 < j < m , is the number of cells whose genotype 
belongs to class j (for instance, having j copies of a specific gene). Then its 
evolution can be described by the system 

d t v = Kv, (1) 

where K is the matrix describing connections between the nodes. For the dy¬ 
namics on the network on Fig. [ 7 J IK can be given by 
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and describe a situation in which the cell in class j divides into two daughter 
cells, one of which has the same genotype as the mother, while the other 
changes its genotype to that of class i with probability kij. In such a case 
( kij)i<i,j<5 cl column stochastic matrix. We note that a particular case of 
this model is the discrete Rotenberg-Rubinov-Lebowitz model where the 
cells are divided in classes according to their maturation velocity. 

Example 1.2. Aristizabal and Glavinovic model of synaptic depression. In pQ 

the authors introduced a heuristic model of synaptic depression. In this model, 
neurotransmitters are localized in three compartments, or pools: the large pool, 
where also their synthesis takes place, the small, intermediate, pool, and the 
immediately available pool, from which they are released during stimulus. The 
key assumption in |T| was that the dynamics of the densities U = (U\, U2, U 3 ) 
of vesicles with neurotransmitters in the pools is analogous to that of voltages 
across the capacitors in the electric circuit. For the electric circuit, E is 
the electromotive source, CiS are capacities, while RiS are the resistances. 
Biologically, E represents synthesis, C^s are the capacities to store vesicles 
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and are interpreted as the pools ’ replenishment rates. This results in the 
following system of ODEs for U, 


dtU = KU + F, 


(3) 


where 
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1.2.2. Macro and micro models. It is worthwhile to reflect on the models dis¬ 
cussed in Examples 11.11 and 11.21 Both represent systems that operate and are 
observable at the macro-scale. This is the scale of our everyday experience 
as, in principle, we can measure the total number of cells with a particular 
genetic characteristics or the number of vesicles in a particular state. Such 
systems, which have characteristics observable at the macro-scale, are called 
macro-systems. The models introduced above only give heuristic relations 
between macro-features, or observables, of such systems, ignoring any under¬ 
pinning dynamics influencing these relations. Models of this kind, which here 
we refer to as macro-models , are relatively easy to use but often are too crude. 
Following the ideas presented in Ref. m, we observe that usually the nodes 
themselves often have internal structure with a complex dynamics, and the 
state of the system in these nodes is the resultant of this dynamics. More¬ 
over, the nodes are interlinked and the rules of connection, given above by 
the coefficients of the weighted adjacency matrix K, may themselves be deter¬ 
mined by dynamical processes taking place along the edges. These scenarios 
can occur in various configurations and the models enhanced by considering 
dynamics at the nodes and/or along the connecting pathways, here called 
micro-models, typically provide a better insight into the dynamics of com¬ 
plex processes. At the same time, in such micro-models the added dynamical 
processes often act at different time, or size, scales and thus the micro-model 
actually becomes a multi-scale model. We observe that the micro-model must 
be asymptotically consistent with the macro-model; that is, the macroscopic 
features of the micro-model should be approximately the same as the fea¬ 
tures provided by solving the macro model. In other words, the observables 
determined by the solutions of the macro-model must be recoverable from 
the multi-scale micro-model by taking its regular or singular limit, and the 
coefficients of the macro-model should be fully determined by the processes 
occurring at the micro scale. We note that such a point of view, though in a 
restricted setting, has been already mathematically explored in CUE]. 

It is, however, important to observe that, in general, building a micro¬ 
model on the basis a given macro-model may lead to a micro-model whose 
singular limit is completely different from the original macro-model. An ex¬ 
ample of such a micro-model is offered by the McKendrick model with ge¬ 
ographical structure and fast migrations between the patches, described by 
a matrix K. It follows HUGO® that the aggregated dynamics in such a 
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Figure 3. Hypothetical micro-model related to a network 
in which the population in the node v q evolves according to 
the McKendrick model influenced by populations migrating 
through the edges e*, and e;, modelled by a functional fj’ 1 . 
On the other hand, the exchange between the nodes occurs 
through diffusion along the edges with fluxes at the end¬ 
points given by some functionals of the population densities 
w p and w q in the nodes and the density Uk on the edge. 


model is given by a scalar McKendrick equation with averaged birth and 
death coefficients, as explained in details in Example 14.31 but not by dU)- 


2. Main results 

Let us consider model (CD, describing the evolution of subpopulations concen¬ 
trated in m locations, with no internal dynamics, which influence each other 
according to the pattern of connections described by the entries of a ma¬ 
trix K. Our aim is to build a micro-model by including internal dynamics in 
the subpopulations, which would have asymptotically the same macroscopic 
characteristics as ©• In this paper we accomplish this by allowing for the 
subpopulations to evolve along the edges of some (hyper) graph according 
to either diffusion or transport operators and subject to specific interface 
conditions at the endpoints of the edges. 

Before we formulate the main results, we have to introduce some basic 
notation. We consider problems in M m and the boldface characters will usu¬ 
ally denote vectors in R m , e.g. u = (m, , u m ). We denote M. = {1,..., m}. 
Further, for any Banach space X , we will use the notation X = X x ... x X, 

7TL tiTTIGS 

e.g. for X = Li(I), I = [0,1] we denote Li(7) = Li(I) x ... x 

V -V-—' 

m times 

The main contribution of the paper consists of the following two results. 
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Diffusion along the edges. Let us consider the following initial-boundary 
problem 


d t u e (x, t) = -d xx n. e (x,t), (x, t) €]0, l[xR + , 

9 x u e (0,t) = eK oo u e (0, t) + eIK 01 u £ (l, t), t > 0, 

9 K u e (l,t) = eK lo u e (0, t) + eK 11 u £ (l, t), t> 0, 

u e (x,0) = u(x), xG]0,1[, (4) 

where e > 0 and u e H = {00,01,10,11} are real mxm matrices. 

Further, let 


/ Ui t€ (x)dx, . . . / U m ,e(x)dx 

Jo J o’. 


and let v be the solution to © with 

itr — nr io nrOO 


roi 


(5) 


and the initial condition given by ^ /q Ui(x)dx ,... u m (x)dxj . Then, for 
any u £ W \{I) and T > 0, there is Ct, independent of e and u, such that 
K0) - v(t)||x < eC T ||u|| W 2 (/) , (6) 

uniformly on [0,T], where X = Li(J) or X = C(J). 

Transport along the edges. Consider the transport equation 

d t u e (x ,t) = ~d x u e (x,t ), sg]0,1[xI + , 

u c (0,t) = u e (l,t) + eBu e (l,t), t > 0, 

u £ (a;,0) = u(a;), x e]0,1[, (7) 

where e > 0 and B is an arbitrary matrix. If v e is defined by m and 
v is the solution to Q with K = B and the initial condition given by 
^u 1 (x)dx,...f^u m (x)dx S j, then, for any u £ W}(J) and T > 0, there 
is Ct, independent of e and u, such that 


||v e (t) - v(f)|| Ll(/) < eC T ||u|| w i (J) , (8) 

uniformly on [0,T]. 

The results presented here heavily depend on the well-posedness theory 
of problems of the form (£!]) and Q and precise estimates obtained for them. 
This theory has been developed in Ref. [5] . Here we only recall these results 
when needed. 


Remark 2.1. The above results show that the models © and © both are 
correct micro-extensions of the macro-model m in the sense explained in the 
introduction; that is, their macroscopic features (here the total size of each 
subpopulation) evolve approximately in the same way as described by (Qp. 

Remark 2.2. We note that (Up and m fall short of typical results expected 
in asymptotic analysis, where one constructs an approximation of the whole 
solution u e by using appropriate initial, and possibly other, layers. This also 
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allows for proving that the singular limit solution v provides a good approx¬ 
imation to u e outside a narrow transient layer close to t = 0, see Ref. [5] . 
We are able to prove such a result for the diffusion problem. However, for 
the transport problem the standard initial layer is a fast-oscillating function 
of zero mean and thus here we cannot claim that the solution to m approx¬ 
imates u c away from zero. 

Remark 2.3. The systems of the form w and & have originated from the 
modelling of transport and diffusion processes occurring along the edges of a 
physical graph, where the exchange between subpopulations only could occur 
at its vertices, see e.g. Refs 

To explain the relation between the latter models with © and ©, let 
us recall that by a graph we understand the pair Q = (V,E), where, say, V 
is the set of n vertices and E is the set of m edges. Note that an edge here 
is an (unordered) pair of vertices. In transport problems, since transport, in 
contrast to diffusion, has a direction, it is more convenient to impose an 
orientation on each edge and work with a directed graph. A directed graph, or 
digraph, is the pair G = (V(G), E(G)) = ({vi,..., v n }, {ei, ..., e m }), where, 
in contrast to Q, the edges here are defined as ordered pairs of vertices. An 
important role is played by the line digraph L{G) of G. To recall L(G) = 
(V(L(G)),E(L(G))) = (E(G),E(L(G))), where 

E(L(G)) = { uv ; u, v G E(G), the head of u coincides with the tail ofv}. 

An analogous definition can be made for the line graph L(Q) of the graph Q. 
To avoid considerations related to the geometry of the edges, see e.g. Ref. eh 
S ection 2.2.1], from the beginning we assume that each edge is identified with 
[ 0 , 1 ]- 

Then, under some technical conditions ®|2S|, the transport of a sub¬ 
stance along the edges of G with Kirchhoff's type interface conditions at the 
vertices can be cast in the form of © with some matrix A describing the 
boundary condition (in ©, A = I + eBJ. However, in such a case, A must be 
the weighted (transposed) adjacency matrix of L(G), see Ref. [7j. 

Similarly, the diffusion along the edges of some graph Q with Robin type 
interface conditions da can be written in the form © but then the matrices 
K w , uj = {00,01,10,11}, must satisfy certain conditions which allow for the 
construction of the (weighted) adjacency matrix of the line graph of Q, see 
Ref. 0. A more detailed presentation of this case is given in Example 1 3.21 

It is easy to see that not every matrix A is a weighted adjacency matrix 
of a line graph. At the same time, there are legitimate models of the form © 
with arbitrary (nonnegative) matrix in the boundary conditions. An example 
is rendered e.g. by the discrete Rotenberg-Rubinov-Lebowitz model [35j . where 
the subpopulations can communicate without any physical connection between 
them. 

We emphasize that the results presented in this paper are valid for arbi¬ 
trary matrices in the boundary conditions of 0 and ©• However, in many 
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cases it is important to determine whether there is a physical graph structure 
behind the model. This question was addressed in Ref. |?J. 

We note that another approach would be to consider hypergraphs as any 
graph is the line graph of a hypergraph ,12], but we will not pursue this line 
of research here. 


3. Diffusion problems 

We shall consider 0 in both C (I) and Li(/). For technical reasons, we 
also shall need the space W((J). The norms in these spaces will be denoted, 
respectively, by || • ||oo, || ■ ||o, || ■ ||i but, if it does not lead to any misunder¬ 
standing, we use X to denote any of these spaces and || • || or || ■ ||x to denote 
the norm in X. Further, by ||| • ||| with appropriate subscript we denote the 
operator norm in respective X. 

To keep the notation in line with Ref. 0 (and with Ref. [ 22 ], on which 
the well-posedness results are based), we introduce operator L by the formula 


X 9 u -> Lu = (7 0 5 x u,7i5 x u) e Y = T x R m , (9) 


where 7 i,i = 0,1, is the trace operator at x = i (taking the value at x = i if 
X = C(/)). The domains of L are D(L ) = C 1 (J) if X = C(/) and D(L) = 
W i(I) in the other cases. Then we define the operator 

X3u-> <Fu = K( 7o u, 7 iu) = ( J^o Ri! ) ( ^” ) e x 


Let A denote the differential expression Au := d xx u. Then we define the 
operators AJ|, a = 00 ,0,1, by the restriction of A to the domains 


D(A%) = {ueC 2 (/) ; iu=$u}, 

D{ A°) = {ueW 2 (l);Lu = $u}, 

D{Al) = (u G W?(/); Lu = $u}, 


respectively. If the base space is clear from the context, we shall drop the 
superscript from the notation. Further, let us denote 


$*u = K*(yoU, 7 iu) 


k oot ^ k iot 

-K 01t K 11t 


7oU \ 
7iu ) ’ 


( 10 ) 


where K T denotes the transpose of K. Clearly 


($*)* = ( 11 ) 

In general, if A is the generator of a semigroup, we denote by {e tA } t >o the 
semigroup generated by A. 
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3.1. Basic estimates 

Consider the resolvent equation for (j4j with e = 1, 

Au — <9 xa; u = f, x £]0,1[. (12) 

Its solution u is given by 

u(ar) = C 1 e-^ + C 2 e' iX + U M (x) (13) 

where 0 A = g? = \\\e ie with Sfyt > 0, 

l 

u ^) = ^/ e~^~Ms)ds. (14) 

0 

and Ci and C 2 are determined by the boundary conditions. Further, denote 
S Q ,e„ = {A = \\\e l6 £ C; |A| > a, |0| < 0 Q < n} 
and So,g 0 AI#q. 

We note that Ao corresponds to the Neumann boundary conditions in 
one of the spaces X = C(J), Li(J), W((J). Standard calculations [20ll8] give 

\\R(\, Ao)/||x < m 2|l/l ' x /9 , A e £* 0 (15) 

| A | cos uq f 2 

for any 9q < 7r. Hence, in particular, Ao generates an analytic semigroup in 

X. 

With regards to J3J in X = C(J), Li(J), W)(I), we have 

Theorem 3.1. [5] Then there is a > 0 such that for any 9 q < tt 

||fl(A, A*)|| x < 2||1?(A, A q )||x, A £ E a ,* 0 . (16) 

Hence, A$ generates an analytic semigroup in X. Moreover, {e tA# }t>o is 
positive if and only if — K 00 , K 11 are nonnegative off-diagonal and — K 01 ,K 10 
are nonnegative. 

For further use we need some intermediate results from the proof. First, 
since in X = C(I),'W\(I), the operator >1) in the boundary condition Lu = 
<Fu is compact (and thus also bounded), Theorem 13.11 is a straightforward 
consequence of Theorem 2.4 of Ref. [22] . More precisely, it follows that 

1?(A,A $ ) = (/-L a $)- 1 1?(A,Ao), (17) 

where L\ = (T|xer(A-A)) _1 : A € p(A 0 ) and the proof relies on the fact that 
if $ is compact and Ao generates an analytic semigroup, then there is a > 0 
such that for A £ S Qj e 0 we have |||I/ A $|||x <1/2 which, by the Neumann 
series, gives (HU). 

The situation with the generation in X = Li(/) is more complicated as 
$ is not bounded, and so cannot be compact, on X. A more involved argument 
shows that the resolvent R( A, A$) exists in the same sector as R( A, A$) and, 
moreover, (i?(A, A|f, ))*f = i?(A, A$)f whenever / £ Li(J) (the adjoint of 
R( A, A|°,) acts in the space of complex vector measures and here we identify 
Li(J) with the subspace of (densities of) absolutely continuous measures). 
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Then, since the norm of an absolutely continuous measure equals the L\ 
norm of its density m , we arrive at 

||-R(A, A^)|| Ll( / } = ||.R(A, Afi)||c(/)- (18) 

Being densely defined, A^, generates an analytic semigroup on Li(J). 


3.2. A lifting theorem 

In asymptotic analysis we need results on solvability of the inhomogeneous 
problem 


d t u(x, t ) 
d x u(0,t) 
d x u(l,t) 
u(x,0) 


dd xx u(x,t) +f(x,t), {x,t) g]0, l[xR+, 

K oo u(0,f) + K 01 u(l,f) + t > 0, 

K lo u(0,f) + K 11 u(l,t) + <pi(t), t > 0, 

u{x), i£]0,l[, (19) 


where f, cf) Q and (f> 1 are known functions and d > 0. First, we observe that 
v(x) = —x(l — x)((a + (3)x — a) satisfies 

a a ,v(0) = K°°v(0) + K 01 v(l) + a, 

d x v(l) = K lo v(0) + K n v(l) + (3. (20) 


Thus, if v is the function constructed above with a = <f> 0 and (3 = <f> 1 then, 
by substitution U = u — v, problem (fTUl) is reduced to 


d t U(x,t) 


d x U(0,t) 
3*U(l,t) 
U(x, 0) 


dd xx \J(x,t) + f (x,t) + 2d((4> 0 (t) + 4>i(t))x - 2 - 0i(t)) 
~x 3 (d t 4> 0 {t) + d t 4>\{t)) - x 2 (2 d t <j> 0 (t) + d t (f>i(t)) + xd t (j) 0 (t), 
K oo U(0,f) + K ol U(M) 

K lo U(0,t) + K n U(l,t), t > 0, 

u(x) + x(l - x)((</> 0 (0) + 0i(0))x - 0 o (O)), (21) 


which is classically solvable, (33], on ]0,T[ for some T > 0 provided e.g. the 
inhomogeneity is locally Holder continuous on ]0,T[. 


3.3. Basic results on solvability of ([4]) 

Let us return to 0 . Let e > 0. We denote by A e ,$ the operator given by the 
restriction of the expression u — > A e u = e~ 1 d xx u to the domain D{ A fi $) = 
(u G X; d xx u G X, Lu = e<f>u}, where X = C(/) or X = Li(/). It is 
easy to see, by m and (ED, that there exists a positive analytic semigroup 
{e tAe , <|, } t >o in X that solves 0 . We show uniform boundeness of {e tAe ’ <I ’} t >o 
with respect to e. 

Lemma 3.1. For any X = C(J), W](7), Li(J), there are constants M^ and 
0!$, independent of e, such that 

||e tA ^||£ ( x) <M^\ 


t > 0. 


( 22 ) 
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Proof. First we observe that since the resolvent equation Au — e~ 1 d xx n = f 
can be written as eAu — d xx u = ef and the Neumann boundary condi¬ 
tion does not involve e, estimate m is uniform in e; that is, for X = 

C(/), L^/), W}(/), 

||i?,(A,e- 1 A 0 )f||x < A G E, 0 . (23) 

| A | cos u o / 2 

Then, as indicated below o. the generation result in X = C (/) and in X = 
W((J) follows due to possibility of establishing the estimate |||La3>|||x < 1/2 
in Yi a j 0 . Here we have to show that for Q such an estimate is valid uniformly 
in e. 

In (j4j) we deal with the family of operators {e _1 A} £> o, hence we have 
to introduce the family La, £ = {L\ Ker ^x_ e -i^)~ 1 . Since Ker(X — e _1 A) = 
I\er(e A — A), we have La, e = L £ a,i = L e \, where the latter refers to the 
case when the diffusion coefficients equals 1, see Subsection 13.11 We see that 
I\er (A — e" 1 A) consists of functions given by u(*) = Cie MX + C 2 e _/ “, where 
fi 2 = eX. Hence Lu = (/xCi — /rC 2 ,/rCie M — /iC 2 e _M ). Thus, for a given 
(di,d 2 ) G R m x R m we have 


L £A (di,d 2 ) 


d 2 — die ^ 
/r(e^ — e - ^) 




+ 


die^-d 2 ~_ ux 
/i(e^ — e _A1 ) 


Then, for X = Li (/) we have 


||L £ A(di,d 2 )|| x 

„ e-*^(||d 2 || + ||di||e-^)(e^-l) e -*^(||d 2 [| + ||d 1 ||e^)(l - e~^) 

~ |/i|3fyi(l + |/*|3ty(l- e~ m ») 

||d 2 || + HdAlle-^ ||d 2 |le-^+||d 1 || C 

- H5ftp(l + e-^) + |/i|^i(l + e-»M) - e|A|cos0 o /2 111 ' u 2,]U 

where C is independent of e and we used the estimate 


1 


OO 

V r~ 2n ^ 

00 — gfy, 

gll _ g —n 


n =0 

- c l _ e -2dtv 

n=0 


(24) 


For X = W}(/) we additionally need to estimate d x L e A(di,d 2 ). Since the 
differentiation only introduces the multiplier ±/i in each term, we have 


||9 a; L e A(d 1 ,d 2 )|| Ll (/) < 


C 

^/e|A| cos $o/2 


II (di, d 2 ) || ■ 


Hence, noting that in (fTTl) the operator $ is replaced by e<h and (di,d 2 ) = 
e<Fu with ||(d 1 ,d 2 )|| < eM 0 ||u|| w i(/) if u G W((J), with M 0 independent of 
e, we obtain 


||eL A , £ $u|| w i (J) < 


Mi 

W\ 


ll u llwj(7); 


(25) 


where M\ is a constant independent of e and A G £ a ,e a some fixed a > 0 . 
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For C(/), we observe that the norm in W}(J) is stronger that that in 
C(/) and also ||(di,d 2 )|| < eM'||u|| C (/) if u S C(/), with M' independent of 
e. Hence 


||eiA,e ( f , u||c(/) < C' , ||eAA,e^u|| w i( J ) < 


M 2 


ll u llc( 7 ), 


(26) 


where M 2 is independent of e and ||A|| and C' is the embedding constant. 
Finally, the result for Li(J) follows from (fl8l) and (l26l) . Thus (l22l) holds. 

□ 


3.4. Asymptotic state lumping 

3.4.1. Formal expansion. We follow the standard asymptotic analysis ap¬ 
proach, PJ. First, we find the hydrodynamic space V of A £i $; that is, the 
null space of u —» d xx u on L>(Ao : o) = D( Ao). This amounts to finding the 
solution of m uncoupled Neumann problems 


(Jxx'Ui — 0, d x Ui(0') — d x Uj(l) — 0. 


(27) 


Hence, V = span{ej}j 62 vi, where e, are versors of K. m . The (formal) spectral 
projection onto V is given by 


r 1 r 1 

Ui(x)dx,... / U m {x)dx 

3 Jo 


v = Vu = 

see ©. Hence we decompose 

U e (x,t) = [PvL e ](t) + [Qu £ ] (x,t) =: V e (t) +W e (x,t), 
where Q = I — V. Further, to shorten notation, we denote 


(28) 


(29) 


^10 


w 

rl 


- K + K - 

10 1 Ttr 11 w 


r 01 


rOl 


r oo 


— jA 11 — K 01 , K° = K 10 — K 00 . 

If u £ is a classical solution to ©, we can apply V to get 

^v e (x,t) — —7^d xx \i e (X)t^ — ~ (*£» t^dX) ..., J* d xx u eni {x^t^dx 

= - {d x u e (l,t) - d x u e (0,t)) 
e 

= (K 10 - K oo )u £ (0,f) + (K 11 - K 01 )u £ (1, t) 

= Kv e (f) + K°_w £ (0, t) + K)_W £ (1, t). (30) 

Similarly, applying Q = (I — V) to ©, we get 

d t vr e (x,t) = d t {u £ {x,t) - v e (a;,i)) = ^d xx (v e (x,t) +w e (a :,t)) - d t v e (a :,t) 

= -d xx , w e (x, t) — Kv e (i) — K°_w e (0, t) — K)_w £ (1, t). (31) 

e 
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The boundary conditions for w e at x = 1 are found from 

d x w e (x,t)\ x=1 = d x u e (x,t)\ x=1 = eK n u e (l,f) + dK 10 u e (0,i) 

= eK+v e (i) + eK n w e (l, t) + eIK lo w £ (0, t). (32) 

Similarly, at x = 0 we have 

<9 x w e (:r,i)| x=0 = eK° v e (i) + eK 01 w e (l,t) + eK oo w e (0,t). (33) 

Now, expanding w e = w 0 + ewi + 0(e 2 ), and substituting it into (13T1) (1531) . 
we find, for e°, 

9xxW 0 = 0 (9 x w 0 (l, t) = <9 x w o (0, t) = 0. (34) 

From this equation we conclude that Wo is constant. Since, by (EH1) . w e and 
thus all terms of its expansion annihilate constants, we conclude that wo = 0. 
Hence the terms containing w e in (1301) are of higher order in e and we can 
consider the following limit equation for v e 

9 t v = Kv, v(0)='Pu. (35) 

Next, at the e level, we have 

<9 xx wi = K v(t), 

9 x wi(l 5 x wi(0,t) = K° v(t), (36) 

where v is given by (1351) . Integrating, using the boundary conditions and the 
definition of W, we find 

Wi(a:,f) = -x 2 Kv(t) + a’K ( |_v(t) + d (t) (37) 

with 

d (t) = - v(t). (38) 

Clearly, the pair (v, ewi) cannot be an approximation of u e = (v e , w e ) unless 
w e | t= o = 0(e). To improve the approximation we introduce the initial layer 
by rescaling time as r = t/e and expand both v(x,t) and w (x,t) in e : 

v(x, t) = v 0 (x, t) + O(e), w(x, r) = w 0 (x, r) + 0(e). 

Standard calculations give Vo = 0 and, using TO, at the e° level, we get 
<9 r w 0 (x,t) = d xx w 0 (x,t), d x w 0 (l,r) = d x w 0 (0,r) = 0, 

w 0 (x, 0) = w(x) = u(x) — Va. (39) 

Since fg xv 0 (x,0)dx = 0, the solution can be written as the Fourier series 

OO 

w 0 (x, r) = e A ” r a n cos(mrx), (40) 

n= 1 

where A n = — (nn) 2 < 0, n = 1,..., and 

1 

a n = 2 J w(x) cos(n7rx)dx. (41) 

o 
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Now, we can formulate the main theorem of this section. 

Theorem 3.2. Let X = C(/),Li(/). Let u e (t) = e tAe ’*u with li G W \{I) be 
the solution of © in X, v be the solution to 13,51) and wo be the solution to 
1,9.91) . Then, for any 0 < T < oo, there is C = C(T, K 00 , K 01 , K 10 , K 11 ) such 
that 

IIu e (t) - v(i) - w 0 (t/e)||x < eC||u|| W 2 (J) (42) 

uniformly on [0, T]. Furthermore, if u G X, then we have 

lirn ||u e (f) - v(t) - w 0 (t/e)l|x = 0 (43) 

£-*■0 

uniformly on [0, T]. 

3.4.2. Proof of Theorem 13.21 The first problem we encounter is that since 
w does not satisfy the boundary conditions of and thus it is not in the 
domain of the generator, w 0 is not a classical solution to (EiiTll up to t = 0. 
To remedy this, we need the following lemma. 

Lemma 3.1. Let u G W \{I) and X = Li(J) or X = C(/). Then for any 
5 > 0 there is 115 G W 1(1) satisfying 115 ( 0 ) = 115 ( 1 ) = 0 and such that 


llu-usllx < Ci(S||u|| W 2 (J) , (44) 

ll^UiUx < C 2 ||u|| W 2 (7) , (45) 

II^u^Hl^/) < C3||u|| W 2 (J) , (46) 

OO 

l a n,<51 < Clllullw^j), (47) 

n= 1 


where Ci,i = 1,2,3 are independent of 5 and a are coefficients ©P for 

u 5. 


Proof. Let u G W 2 (J). Then d x u absolutely continuous and, for any x G I, 
we have |u(x)| < C'||u|| W 2( / ) and |<9 x ii(x)| < C"||u|| W 2 ( 7 ). For 0 < S < 1/2, 
we define 


( u O) 

for x G [<5,1 — A], 


U 5 (x) - { X J ) {x 2 £ 2 ) + u(<5) 

for x G [0, <5[, 

(48) 

[_ —-(<5 2 — (1 — x) 2 ) + u(l - 5) for x G]1 — 5, 1]. 


We see that 



r <9a-u(x) 

for x G [<5, 1 — <5], 


d x us(x) — < X 5 () x 

for x G [0, 5[, 

(49) 

{ (1 x 

) for xG] 1 — 5,1] 


and, since d x 115 is continuous at x = <5,1 

-5, 


( d xx u(x) 

for x G [5,1 — 5], 


d xx \i 5 (x) - < x s ( ' ) 

for x G [0, 5[, 

(50) 

1 — 9 x u(l — S) 

' S 

for x G] 1 — 5,1] 
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hence 115 £ W \{I). Now, we have 


OOO 

l u - u <5||l!(/) < J |u(x)|dx + J \x 2 - 6 2 \dx + u(5) j dx 


+ 


<9zu(l - 6 ) 


25 

1 -s 


1 1 
J |(1 — x) 2 — 5 2 \dx + u(l — 5) J dx 

-S 1 — 5 


and 


|U U <5 1| C(I) < sup |u(x) - u(<5)| + sup 
a?E[0,<5] ccG[0,<5] 


dxU ^.(x a - 5 2 ) 


25 


+ 

sup 


xG[l—<5,1] 

+ 

sup 


a:G[l—<5,1] 


d x u(l - 6) 


25 


((1 - x ?~ 5 2 ) 


< C*! <5 11 u 11 14^2 (, 


where we used |u(x) - u(y)| < sup ?e/ |3 x u(£)||z - y\ < ||u|| W 2 (/) |a; - y |, 
x,y £ I. This gives (l44l) . Analogous calculations, using (l49l) and (l50l) give 
(TI§1) and m- 

Finally, integrating twice by parts and using the boundary conditions 
yields 


1 

2 


^n,d 


- 5 — ~d x u s (x) cos nnx\ x x= l 
'K A n. 


1 

7 r 2 n 2 


1 

J dxxUg(x) cos mrxdx. 
0 


Hence 


and thus JED holds. 


|a„,a| < 


\d xx u s (x)\dx 


□ 


Further, we observe that if 115 is such an approximation of u, then 
w = u — Vvl, where V is defined by (1251) . can be approximated in the same 
way by wj = 14 — Vu g. Indeed, since Vu g is a constant, w 5 is in Wf (/), 
d x wj = so that the boundary conditions are satisfied and 

1 

w - wj = u - u s J W*) — ug(x))dx (51) 

0 

so that wj approximates w in each X: 


||w - w 5 ||x < I|u - u,5 1| x + I|u - u 5 || Ll(7) < 2C'i<5||u|| W 2 (j) . 

In what follows, we shall use the initial layer W 0.5 which is the solution to 
(1551) with the initial condition wg defined as W 5 = U 5 — Vug. We approximate 
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(v £ ,w e ) by (v, ewi + w 0 i j) and find that the error e v = v £ — v, e w = 
w £ — ewi — wo ,5 formally satisfies, by ©, m and and the equations 


d t e v (t) 


dt.Q-W (*£, t') 


e v ( 0 ) 

e w 0 ) 


Ke v (t) + K)_(e w (l, t) + ewi(l,f) + w Q , 5 (l,r)) 
+K°_(e w (0,t) + ewi(0,t) +w Oi 5 ( 0 ,r)), 

ECe v (t) 

-K]_(e w (l, t) + ewi(l, t) + w 0i(5 (l, r)) 

—K° (e w (0, t) +ewi(0,t) +w o , 5 (0,r)) - ed t wi, 
v e (aj,0) - v(a;,0) = 0, 

-ewi(x, 0) + w(a;) — w^(0). (52) 


For the boundary conditions, using © and we have 
3 x e w (l,t) = d x w £ (l,t) - ed x wi(l, t) - d x -w 0 „s(^,r) 
= e (K^e v (t) + K n e w (l, t) + K lo e w (0, t)) 
+ e(K 11 wo,5(l, r) + K lo w o , 5 ( 0 , r)) 

+ e 2 (K 11 wi(l,r) + K 10 wi(0, r)). 


Similarly, using and again we obtain 

3 x e w (0,f) = 3 x w £ (0 ,t) - e^wi(0 ,t) - 3 x w 0 ,s(x,r) 
= e (K° e v (t) + K 01 e w (l, t) + K oo e w (0, *)) 
+ e(K 01 w 0 , 5 (l, r) + K oo Wo,i(0, r)) 

+ e 2 (K 01 wi(l,r) + K oo wi(0, r)). 


To proceed, we observe that, since the semigroup {e tA ‘^} t >o is analytic, for 
a given u £ X the function t —» u £ (f) = e* Ae ’ # u is a classical solution of the 
problem in the sense that for any t > 0 it is infinitely (classically) differen¬ 
tiable and the system of differential equations and the boundary conditions 
are pointwise satisfied. Similarly, all other terms of the asymptotic expansion 
are sufficiently regular for all terms in the above equations to be well defined. 
Hence, denoting E = e v + e w , adding the equations for e v and e w and using 
the regularity of each term, we see that E is the classical solution of the 
problem 


5 t E 

= -d xx E - e9 t wi, 
e 

H 

H 

II 

o 

= eK 01 E| x= i + eK 00 E 

dx^\x=i 

= eK n E| x=1 + eK 10 E 

E(0) 

= — ewi(0) + w — wj, 


x—0 + £(^0,5 + e ^o)> 
x=0 + e(’® r i,<5 + e ^i)> 


(53) 


= K 01 w 0 ,5| x= i + K 00 w 0i 5| x=0 + e(K 01 wi| x= i + K 00 wi| x = 0 ), 
= K 11 w 0 ,i| x= i + K 10 w 0 ,5| x=0 + e(K U wi| x= i + K 10 wi| x = 0 ). 


where 

* 0,5 + £*0 
*i,5 + £*i 
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Next, by m, we see that the substitution 
F (x,t) = E(x,t) + ex(l-x)(('&o,s( T ) + i &iAT) + e (®o(t) + &i{t)))x 

-*s(T)-e* 0 (t)) 

reduces m to 

9 t F = Ia ra F + eG + G 5 , 

6 

a,F | x=0 = eK 01 F | x=1 + eK 00 F| x=0 , 
d x F \ x=1 = eK n F | x=1 + eK 10 F| x=0 , 

F(a;,0) = w(x) - wg(x) + e(-w 1 (x, 0) + H 5 (a:, 0)), (54) 

where, since wi( 0 , 0 ) = — v and wi(l, 0 ) = ) v, 

we have 

G (x,t) = -dtWi(x,t) + 2{&o{t) + &i(t))x-2&o(t)-&i(t), 

G s(x,t) = 9 r (-s 3 (f 0 ,iW + ^i,«(r)) 

- x 2 (2'F 0> s{t) + ^fi,s(r)) + x& 0 ,5 (t)) , 

H 5 (a;,0) = x(l - x) ((k+w 5 (1) + K)J_w 5 (0) 

+ lv(K 01 +K 11 ) (K} + + 2K°) - (K 00 + K 10 ) (K° + 2K^)) x 

- (K 01 W5(1) + K OO W5(0) 

(K 01 (K^ + 2K° + ) - K 00 (K° + 2K^)) v)) . 

We observe that is well defined as for r > 0 


w 0 ,5(0, t) = ^e A " T a n<s , w 0 ,5(l, r) = ^(-l) n e A " r a„, 5 

n =1 n =1 

and, by Lemma |3. II both series are uniformly convergent on [0,T],T < 00 , 
and thus define continuous functions on [0, T], Thus the values w Oi a(0, 0), wo,a(l, r) 
are well defined and equal to wa( 0 ) and wa(l), respectively. 

Similarly, we observe that for r > 0 


dr^o^r) = K 01 ]T(-1)" A„e A " T a M + K 00 A„e A " T a n , 5 

n =1 n =1 

and, by (|T71) . 

°P °»OC X 

J \d T ^o t s(r)\dt < C J ^ 7 me _7r ^ \^ n ,s\ dt I 8 ".*! - eC 4||u|| W 2 (/) . 


0 0 

In the same way 


\d T ^i,s(T)\dt < eG 4 |]u|| W 2 (J) . 
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Then, using the mild formulation of the solution to (1541) . we obtain on [0,T] 
||F(i)|| x < ||e tAe -*(w -w { + eH 5 )|| x 

t 

+ || J e^ A ^(eG(s) + G(s/e))ds\\ x 
0 

||G(s)|| x ds+ J ||G«(a/e)|| X dsj 

< (||w-W 5 ||x+ eC"||u|| W 2 (/) j 


< M$e“* T l\\w-w s \\x+e\\H s \\x+eJ 


for some constant C'. Using the fact that C' is independent of 8 and the 
smallness of the auxiliary terms, we see that for any T < oo there is C such 
that 

IIu e (t) - v(t) - w 0 (f/e)|| x = eC||u|| W 2 (J) 

uniformly on [0, T], where u £ Wf(I) is arbitrary. Using the density of 
W 1(1) in X, uniform boundedness of { e tA e,®} t > 0 0 n [0,T] and corollary 
to the Banach-Steinhaus theorem (3 — e lemma), we also obtain (1431) for any 
u G X. 


Example 3.2. It is interesting to consider to which extent the properties of 
the original system are inherited by the limit one. For this we return to the 
neurotransmitter model of Aristazabal and Glavinovic, described in Example 
m This macro-model was enhanced to a micro-model along the lines de¬ 
scribed in the Introduction in Refs [HI ITS] . where the authors corrected the 
first attempt of such a construction developed in Ref. HU. In the micro-model 
it is assumed that the pools with vesicles occupy some physical space and the 
vesicles in each pool move according to a diffusion process and migrate be¬ 
tween the pools by crossing semipermeable membranes separating the pools. 
In Ref. 1 14] , the authors interpreted the pools as the edges of a finite graph 
connected by the nodes at which the exchange of vesicles takes place. In this 
particular model the graph is linear, however, in Ref. [15] the author general¬ 
ized the model so that it was posed on a finite graph without loops Q = (V, £) 
with, say n vertices and m edges. On each edge there is a substance with den¬ 
sity Uj,j £ A4, which diffuses along this edge and can also enter the adjacent 
edges across the vertices that join them according to a version of the Fick law. 
To write down its analytical form, first we note that, since diffusion does not 
have a preferred direction, we can assign the tail, or the left endpoint, (that 
is, 0) and the head, or the right endpoint (that is, 1) to the endpoints of the 
edge in an arbitrary way. Let lj and rj be the rates at which the substance 
leaves through, respectively, the left and the right endpoints and Ijk and 
rjk be the rates of at which it subsequently enters the edge . Then the Fick 
laws at, respectively, the head and the tail of e, are given by 

- d x Ui( 1) = r-jUj(l) - ^2 r ijUj(v), d x Ui( 0) = ZjUj(O) - ^ kjUj(v), (55) 
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where we used Uj ( v ) as v may be either the tail or the head of ej. If there are 
no edges incident to ej, then respective coefficients rij or Ijj are equal to 0. 

It is clear that if r^ 7 ^ 0, then lij = 0 and if hj 7 ^ 0, then rij = 0 and 
thus we can define an m x m matrix A by setting atj = 1 if either rij = 1 
or lij = 1 and zero otherwise. Then A is the adjacency matrix of the line 
graph L{Q) of Q. Such a matrix is, however, not easy to use and we see that 
introducing, for any i,j £ A4, 

= -kj if v = 0, k°j - -hj if v = 1, fc£° = h, 

kjj = r^ ift> = 0, klj = r^ if v = 1, k\l = —r it (56) 

where v is either the tail or the head of the edge under consideration, the 
boundary conditions can be written as in {^]). In Ref. |15l the author consid¬ 
ered the expected position of the particle on the graph which led to a Feller 
process and thus the problem was posed in the space of continuous functions. 
On the other hand, in Ref. El the probability density of the process was con¬ 
sidered and this resulted in a Markov process posed in a space of integrable 
functions. The corresponding diffusion problems are adjoint to each other, 
as in m and m- Certainly, this interpretation requires in both cases the 
solution to be a scalar function and thus we identify C (I) (resp. Li (I)) with 
the space (7(1) (resp. Li(T)), where I = [0i, li]U.. . U[O m , l m ]; that is, instead 
of considering a vector function on I we consider a scalar function on a dis¬ 
connected compact space composed of m disjoint closed intervals, see 0 . In 
particular, each edge ej,j £M, is identified with the closed interval [Oj , 1 j]. 
Then A$ will be changed to A$ which is the restriction of 

(■ Au)(x ) = Y X[Oj,i j](x)cxjd xx u(x) 
jeM 

to the space ofC 2 (rep. Wf) functions on Int I =]0i, li[U ... U]0 m , l m [ which, 
in the former case, are C 1 on I and satisfy 

mm mm 

d x u(0j) = y k jk u (°k)+Y k % u 0-k)’ d Mh) = Y k j° u (°k)+Y k )l u ( i k)- 

k =1 k— 1 k— 1 k =1 


Denote by {e tAe ’*}t >0 such a realization of {e* Ae ’ <& }t >0 in X = (7(1) or 
X = Li(\). Recalling the relation between the model of Refs mm and our 
formulation, let <f> correspond to the matrices K“, to £ {00,01,10,11}, accord¬ 
ing to \56\) . Then, by Ref. [S], both {e tAc ’*}t >0 and its adjoint {e tAe ^'}t> 0 , 
see m, are positive, in (7(1) and L i(I), respectively, for any e. Let us con¬ 
centrate on the process in Ai(I)- It is clear that it is a Markov semigroup if 
and only if 


0 



2 = 1 


-u( 0 i)Y( k ' 

J=i 


+ k ij) + u(li) Y;( k ij + k 


uOll 


1=1 
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As the end-point values are arbitrary, {e 4A '.' 1 '* } t > 0 is Markov if and only if 

m m 

J2(k%+k™) = 0, B fc «+*& 1 ) = 0. i&M ’ ( 57 ) 

i=i j =i 

which, according to \56\) gives U = ]T) ■ hj an d r i = X ]j r ij > as stated in 

Ref. [T5] . 

If we consider this formulation in the approximation \35\) 
d t v = Kv, v(0) = Vii, 

then 

K = -K 10T - K oot + K 11t + K 01t . 

By [M), K is positive off-diagonal and, by m it is a Kolmogorov matrix and 
thus i f 35]) represents a continuous time Markov chain in which the original 
distributed states are ‘lumped’ together at the vertices. 


4. Transport problems 

4.1. Solvability of the general transport problem 

In the same way as with diffusion, we begin by considering 0 with e = 1 
and denote then T = I + B. The considerations of this section, however, are 
valid in a more general setting [£ . Hence, consider 

u t (i,l) = [Au](i,i), u(0,t) = Tu(l,t), u(a;, 0) = u(x), (58) 

where A is the realization of the expression Au = —9 a u on the domain 
D( A) = {u £ Wj(/); u(0) = Tu(l)}. We emphasize that T is an arbitrary 
(not necessarily nonnegative) matrix. For simplicity we only consider this 
problem in X = Li(/), since in C (I) the domain D( A) is not dense. Then 
we have the following result (5J. 

Theorem 4.1. The operator (A, D(A)) generates a Co-semigroup on Li(J). 
The semigroup is positive if and only if K > 0. 

Denote by {e tAT } t >o this semigroup. Then, [51125], 

[e tAT u](a;) = T"u(n + x — t), —n<x — t<—n+ 1. (59) 

4.2. Asymptotic state lumping in transport problems 

It is easy to see that integrating (1551) and denoting 

v(f) = ui(x, t)dx,... J u m (x,t)dx'j , 

as in (1251) . we obtain the system 

d t v{t) + (I — T)u(l,i) = 0. 
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If one could assume that u slowly varies along each edge; that is, u(x, t) ~ 
constft) = u(l,f), thus giving v(t) = fg u(x,t)dx « u(l,f), then v could be 
approximated by the solution v of the system 

d t v = (T - I)v, 

as in the approximating equation (13511 for the diffusion equation. In the lat¬ 
ter case the above assumption was physically plausible as diffusion tends to 
flatten fluctuations of the density and, with fast diffusion, the density quickly 
becomes homogeneous in space. With transport, the situation is not so ob¬ 
vious as one of the characteristics of transport processes is that the initial 
profile propagates in time with little or no distortion and the choice of a suit¬ 
able scaling is not obvious. In Section f4.2.4l we will show some cases which 
do not yield the expected result but here we will continue with a properly 
rescaled model. 

4.2.1. Properties of solutions to (ff|l. For each e > 0 we define an operator 
(A e , D( A e )) as the realization of the expression A e u = — -d x u on the domain 
D(A e ) := {u £ W)(/); u(0) = u(l) + eBu(l)} and the limit operator 

(A, 15(A)) as Au := — d x u restricted to D(A) := (u £ Wj(/); u(0) = u(l)}. 
Notice that, as for the diffusion, each operator from A e , e > 0, is defined on 
a different domain. By Theorem 14.11 for each e > 0 there exists a semigroup 
{e tAe }t>o- We also have the ‘limit’ semigroup {e tA } t >o- As in the diffusion 
case, we show that {e 4Ae }t>o are bounded uniformly in e. 

Lemma 4.1. For every eo > 0 there is a constant M such that for any e G]0, eo] 
and t £ R + we have 

||e Ae *u||x < Me^ B ^[|u||x- (60) 


Proof. In our case, (Eg becomes 

[e A '*u](x) = (I + eB)”u + x — for — n < x — —t < —n + 1. 

Hence, for n— 1 <t/e<n we have, by changing variable of integration, 


7-n+l 


|e A ' 4 u|| x < ||(I+ eB ) n \\J 


u | n + x —- 


dx\ 


+ ||(I +el)"- 1 1|^ 

\-n +1 

<(l + e||B||n|u|| x . 


U ( 71 — 1 + X- 


dx 


( 61 ) 


As n < | + 1, 

(1 + e ||B||r ||u|| x < ((1 + c ||B||)^OT) l|B||t (l + e ||B||) ||u|| x < Me*" ||u|| x . 


□ 
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4.2.2. Formal expansion. Formal steps are similar to that in Section 13.4.11 
so that we only provide a brief summary of them. First, the hydrodynamic 
space V consists of solutions to 


<9 a .u = 0 , u( 0 ,i) = u(l,t), 

so that V is spanned by {e.;}j e yK, where ej, j G A4, are versors of R m . 
Clearly, the formal adjoint problem is given by 


d x <f> = 0 , 4>{0,t) = 0 ( 1 , t), 

and thus the projection onto V is given as in ©; that is, by 


Vu = 


1 1 

ui(x, t)dx, J U 2 (x, t)dx, ..., J u m {x,t)dx 


(62) 


vo o o 

Projecting the solution of m onto the hydrodynamic and kinetic subspaces 
we have u £ = Pu f + Qu £ = v £ + w £ , where Q = I — V so that w £ G W = 
{w G X; w {x)dx = 0}. Projecting the equations onto V and W, we get 


9 t v e (t) = Bu(0,t) = B(v £ (t) + w £ (0,t)), v £ (0) = Vu 


(63) 


and 


d t w £ (a;,t) 
w £ (0 ,t) 
w e (a:, 0 ) 


— -d x w e (x,t) — B v e (x,t) -Bw e (0,t), 

w e (l,t) + eBw c (l, t) + eB v £ (t), 

u(a;) — Pu. (64) 


Substituting the expansion w e (i r,t) = w o(x,t) + ewi (x,t) + 0(e 2 ) into (IMll . 
first we find that wq = 0. Thus the limit problem for v £ is given by 


d t v(t) = B v(t), v(0) = Vu. (65) 

For Wi (which only is needed for a technical reason) we obtain 

d x wi = — Bv. ( 66 ) 

Thus wi(r,l) = Mv(t)x + H(t), where T~L(t) is arbitrary. Since Wi G W, we 
find wi(i ,t) = Bv(t) (1/2 - x), where v(t) = e Bt Vu. 

As in the diffusion case, the approximation (v £ ,w £ ) ss (v,ewi) does 
not hold close to t = 0 unless u = Vii + O(e). Thus, as usual, we introduce 
the initial layer v(r) = v 0 (t) + 0(e) and w(x,r) = w 0 (a:,r) + O(e), where 
r = t/c. Then standard argument gives Vo^t) = 0, while for the kinetic 
part we get 

d T w 0 (x, t) = -d x w 0 (x,r), w 0 (0, t) = w 0 (l, r), w 0 (a;, 0) = w = u(x) - Vii. 

(67) 

We observe that, denoting by (e" A, ) t >o the semigroup solving (1671) . we can 
write 

[e 7Al w](a;) = w(n + x — e -1 f), — n < x — e~ 1 t < —n + 1 

and, by estimate G2J, 


||e« Al w|| x < ||w|| x . 


(68) 
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4.2.3. The error estimates. Let us define approximation as follows 

(v £ (t), w e (x,t)) = (v(t),ewi(x,i) + w o(x,t)) + (e v (x,t),e w (x,t)). (69) 

As with the diffusion, we encounter the problem that, in general, a given 
initial condition u, even if it belongs to W)(/), it does not belong to D( A e ) 
for all e > 0. Thus the situation is even worse than in the diffusion case, as 
the transport semigroup is not analytic and thus we do not a differentiable 
solution for t > 0, which is essential for the error estimates. Similarly, in 
general, the initial condition in in is not periodic so that the initial layer 
is not differentiable either. Thus, as before, we will work with approximate 
initial conditions. 


Lemma 4.2. Let u £ W](J). 

Qua = vs + Wj £ W](7) that 

Then, for any S £]0,1/2[ there is ua 
satisfy 

= Vu.5 + 

u«(0) 

= ua(l) = 0, 

(70) 

wa(0) 

= wa(l) =0, wa £ W, 

(71) 

ll u <5 ^ u llx 

< C5||u| w i(/), 

(72) 

11 U (5 11X 

< C u w i(/), 

(73) 


for some constant C independent of 5. 


Proof. Let u £ W)(J). Condition (1701) can be achieved by a construction 
similar to that of Lemma o One easily checks that ua defined for 0 < <5 < 
1/2 as 


I 

u(x) 

for 

x £ [d, 1 — <5], 


ua(x) = < 

, f 2 * 

for 

x £ [0, 5) , 

(74) 

1 


for 

x £ [1 - 6,1 ], 



satisfies ua £ W((J). For arbitrary x £ [0,1] we have 

|u<s(aO| < C || u llw}(/) (75) 

and 

||u - u 5 || x < 5 - [u(<5) + u(l - <5)] < C'S ||u|| w i (J) . 

However, for [Qua] we only have [Qua](0) = [Qua](l); that is, the projection 
of ua onto W is periodic but does not satisfy (1711) . Thus, we define 


ua(x) = vs(x) + w 5 (x) 


us{x)dx 



los(x) / u s {x)dx 


, (76) 


f ll i 

\ (l-S)S 

fi 1— X 

{ TT^T5 


x £ [<5,1 — <5], 
x £ [0,<5], 
x £ [1 — (5,1]. 


where 


for 

for 

for 


(77) 
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Since f 0 u>g(x)dx = 1, we have wj £ W. Then we have 


|u-uj|| x < \\Vu - VusWx + ||u - Ui||x + 

1 

< 2 C , ( 5 ||u|| w i (/) + J ||u(a;) - Ui(a:)|| x da; 

0 

1 

+ II1 - wallx J |us(o;)|da: 


1 

1 

/ u (x)dx — u$ 

/ u s(x)dx 

j 

0 

J 

0 


< 3C' , 5||u|| w i(j) + 2(7" <5(1 - <5)||u|| w i(j) < C'i5||u|| w i(_ r p 
Inequality Cl follows by integrating d and using d- □ 


Now we can formulate the main theorem. 


Theorem 4.1. For any T £ ]0, oo[ there exists (7(T, B) such that for any 
(sufficiently small) e > 0 and u £ W((7) the solution u c (t) = e tAe ii of [&!)) 
satisfies 

l|u e (i) - v(f) - w 0 (i/e)|| x < eC(T,B)||u|| w i (/) , (78) 

uniformly on [0, T], where v and wo solve, respectively, I6'5I) and j6'7 [ ). 

Proof. Let us now consider 0 with u £ W)(/) replaced by = (vj,wj) 
constructed as in Lemma 14.21 and consider the corresponding approximation 
of u e 5 = (v e< 5 , w e ,s) given by (vj(t), ewi^x, t)+w 0 s{x, r)). As before, the error 
is defined by e v s(t) = v £ s(t) - vg(t) and e w s(x,t) = w £ s(x,t) — ewq s(x,t) — 
Wo s{x, t). Thanks to the construction of u^, all terms in the error equation are 
differentiable and satisfy the boundary conditions for t £ [ 0 , oo[, in particular, 
we have w O 5 ( 0 ,r) = 0. Hence, by direct substitution, using (1M1) and (l67l) . we 
find that the error is a classical solution of the following problem 


IdtOvg (t) 
dte w g {x, t) 

^w<5 ( 0 ; 

{x , 0 ) 


Be v 5 (t) + Be w a(0, t) + eBwi^O, t), 

— d x e wS (x,t) -B e vS (t) -le W(5 (0,f) - eBw ls {0,t) 
e 

-ed t wi s{x,t), 

e 2 

e w5 (l,t) + eBe v<5 (t) + eBe W(5 (l,f) + —M 2 v 5 (t), 

0 , e w5 (x,0) = -6Wi(i, 0). 


Combining these two equations for the single error E$ = e v ^ + e w ^, we get 


dt'Eg = — d x Es - edtwia, 

e 

e 2 

E,s(0, i) = Ea(l,f) + eBEa(l,i) + — B 2 V 5 (t), E^x, 0) =— 


(79) 


£Wij(s,0). 
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To be able to use the semigroup estimates for (1791) . we have to lift the inho- 
mogeneity in the boundary condition to the equation. To do this, we define 
Fa( 2 :,£) = Fig(x,t) — e 2 (l — x)M 2 vs(t,)/2. Then F 4 satisfies 

d t F s(x,t) = ~^d x F s (x,t) + eJ(x,t), 

F 5 (0 ,t) = F s (l,t) + eBF 4 (l,£), Fj(a:, 0) = eH(z), 

where 

J(ar,£) = ^Bv 5 (£) - |(1 - x)M 3 v s (.t) - Q - x'j B 2 v 5 (£), 

H(x) = - ^ - x^j Bv 5 (x) - J ^ lvg(x). 

Thus, using the Duhamel formula and ([60]l . we obtain 
l|Ei(*)||x < ||F 5 (t)||x + —||B|| 2 e l|B||T ||v 5 || x <— - + ^ ^^ ||vg||x 

eM(e ll B H T — 11 

-“(2 + e||B|| 2 + 2 ||B||)||v 4 || x 

+ jl|B|| 2 e l|B||i ||v 5 || x < Cfe\\vs\\x. 

for some constant C' T independent of e and the initial condition. 

To complete the proof, we take u £ W}(/) and, for any given e > 
0 we construct the approximation with <5 < min{e, 1/2}. Further, let 
v + ewi + w 0 be the approximation constructed by (l65l) . (l66l) and (l67ll with 
v = ?ii, w = Qii. Then, denoting by E the error of this approximation, we 
find 

||E(t)||x < ||e* A 'u - e* B v - ewi(£) - w 0 (r)||x 

< ||e iAe u 5 - e m -vg - ewi 4 (£) - w 04 (t)|| x + ||e* Ae (u - u 4 )|| x 

+ ||e tB (v - v 4 )|| x + e||wi(£) - w w (£)|| x + ||w 0 (r) - w 0(S (r)|| x 

||]g|| ji 

< C're||v 5 || x + Me l|B||T ||v - v 4 || x + ^—||v - v 4 || x 
+ ||w — w 4 || x < eCr||u|| w i(j), 

where in the last line we used <E3>, C3 and 5 < e. Then, noting that 
e||wi(£)|| x < ee T ll B H/2, we find that (1751) is true. □ 


Remark 4.1. We observe that this result is of a different type than Theorem 
\3.2\ In the latter, the initial layer term Wq decays exponentially to 0 as e —> 0 
for any t > 0, see 1 7W - Indeed, let t £ [to,T],to > 0, 0 < e < eo- Then 
can be estimated as 


w 0 (£/e) < Ci(£ 0 ,e 0 ,u)e 


Ci(£ 0l e 0) u) n 2 t _ 2 £ 1 
< e- 7 - -e « < C 2 e 

7 T Z tn £ 


where C 2 is independent oft and e. Hence can be written as 


u ft) - v(t)|| x < eC 3 
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uniformly on [to,T], where C 3 is independent of e and t, but depends on u, 
to and the coefficients of the boundary conditions. In other words, outside an 
0(e) transition zone, the whole solution to the PDE problem on a network 
w can be approximated by the solution of an ODE system 1.9,5)) . 

This is in contrast to m where the term wp does not decay exponen¬ 
tially with e —> 0 but is periodic with period e (in t). Hence, in the transport 
problem 0, the solution u e (f) cannot be approximated by the solution v to 
the Cauchy problem for an ODE, defined by However, using the fact 

that X = V © W, we can project to find that 

||v e (t)-v(t)|| x < eC(T, B)||u|| w i (J) , ||w e (t)-w 0 (V e )llx < eC(T,B)||u|| w i (J) . 

(80) 

Thus, the macroscopic characteristics of the flow on the network; that is, the 
mass on each edge, can be approximated by the solution of ODE \65\) . The 
reminder, however, approximately behaves as a fast oscillating function with 
zero mean. 


Remark 4.2. In the simple case of unit speeds along each edge, the convergence 
ensured by the first estimate of TM ) can be proved directly. Recall that 


[e Act u](a;) = (I + el 


( n + x — -t 


for — n < x - 1 < — n + 1 . 

e 


Hence, for n — 1 < t/e < n we have, as in l6l\) . 


1 1 

'P[e Ae *u] = (I + eB)™ _i yu (z) dz + eB(I + eB ) nl J u ( z ) d z ■ 


We have t = (n — l)e + t', where t' G [0, e]/ that is, n — 1 = t/e + 9(e) with 
0 < 9(e) < 1. Since B is a matrix, using the Dunford functional calculus we 
obtain 

lim(I + el ) 71 " 1 = lim(I + eB) i (I + el) 9(e) = e tB , lim eB(I + eB )"" 1 = 0. 

e —>0 e —>0 e —>0 

Hence 

lim V fe Aet ul = e tB Pu, 

e ->0 1 J 

where e tK V\i is the solution to IU5\) . We observe that this result was obtained 
without assumption that u £ W((/) but also it does not yield the rate of 
convergence ensured in (GH). In fact, using the density of W((/) in X and 
uniform boundedness with respect to e of all involved operators, this result can 
be deduced from \80\) by the 3 —e lemma (a corollary to the Banach-Steinhaus 
theorem), as in 0 - 

4.2.4. Some pitfalls of constructing micro-models. In the first two examples 
we consider the transport model (l58l) on a strongly connected digraph G so 
that the matrix T is irreducible and column stochastic |3) [22]. 

Example 4.1. First we try to mimic the scaling of 0 and consider 

d t u e + = 0, u f (0, t) = eTu e (l, t), u e (a;, 0) = u(x). (81) 

e 
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It is clear that the hydrodynamic space of this problem, spanned by the solu¬ 
tions d x u = 0, u(0,t) = 0, only consists of the zero function and thus 1181(1 
does not offer any interesting limit dynamics. 

Example 4.2. In the next case the boundary condition is changed as follows, 

dtu e + ^-d x \i e = 0, u £ (0, t) = Tu c (l, t), u e (a;,0) = u(x). (82) 

In this case, the hydrodynamic space V consists of solutions to d x u = 0, u(0, t) = 
Tu(l,f). Proceeding as in previous section, we find that the projection of u 

onto V is given by Pu = ^1 • f^ ufajdx'j N, where N is the Perron eigen¬ 
vector of the matrix T and 1 = (1,..., 1) is its left eigenvector; we assumed 
that N is normalised so that 1 ■ N = 1. Then v e = V u e satisfies 

d t v e = ^ Pd x u t = i ((u(l, t) - u(0, £)) • 1) N = ^u(l,t) • (I - T T )1 = 0. 

Hence, the equation for the projection of the solution onto V (the hydrody¬ 
namic part of the solution) is exactly the limit equation and v e = pN. where 
p = fo Ui(x)dx^ is the initial mass. Then it is easy to see that the 

complementary projection (the kinetic part) w coincides with the initial layer 
w. However, the initial layer equation is identical to \82 1) (with e = 1) but 
with the initial condition in the kernel of V. The asymptotic theory of such 
problems has been developed in e.g. Refs. [3j 126} [30] and, in general, it follows 
that w does not have an exponential decay as t ^ 0 (e 0). We obtain the 
(exact) decomposition 

u e (a;, t) = pN + w(x, t/e) 

which somehow resembles Theorem EH However, the hydrodynamic part of 
the solution gives the stationary distribution of the total mass on the network 
among the edges. Hence, it provides neither simplification of the original prob¬ 
lem nor any interesting approximate dynamics. 

Example 4.3. Structured McKendrick model (?) Tj [o( ;(T|. Consider a popula¬ 
tion divided into m groups with respect to some attribute i and assume that 
the individuals can move between them. These groups could refer e.g. to ge¬ 
ographical patches (and then the change of the group would be due physical 
migration), or to the number of particular genes (whereupon the change can 
occur due to mutation). If the vector n (t) = (jii (£),..., n m (t)) £ R m gives 
the numbers of individuals in groups 1 ,... ,m at time t, then 

d t n = Kn, n(0) = n, (83) 

where K = {/cr/}i<i,j< m is a Kolmogorov transition matrix (of a time-continuous 
process); that is, nonnegative off-diagonal and the columns sum up to 0. As¬ 
sume now that we also want to include the demographical processes in each 
patch and that migrations occur at a much faster rate than the demographic 
processes. Let us denote by the population density in patch i and by a the 
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age of individuals. Then in each patch the density ni tE (a,t) should satisfy the 
McKendrick equation and hence we arrive at the system 

OO 

d t n e = An e + Mn f H—In f , n e (0,i) = j M(a)n e (a,t)da, n(a, 0) = n, 

o 

(84) 

where A = diag {— d a , ..., — d a }, B(a) = diag{/3j(a)}i< :; '< m , describes the age 
and patch specific fertility rates, M(a) = diag{— Pj(a)}i<ij< m , where Pj(a) 
is the age specific death rate in patch j and 1/e is the ratio of the reference 
times of demographic and migration processes. 

Note that here we have a simplified version of the model from Fig. [21 
where the exchange between the nodes occurs instantaneously according to the 
migration rates kij, 1 <i,j < m,i ^ j. 

Biological heuristics suggests that no group structure should persist for 
very large intergroup transition rates; that is, for e —> 0. The precise re¬ 
sult depends on the structure of the network of connections described by the 
matrix K. In particular, if the corresponding network is strongly connected; 
that is, if IK is irreducible, then A = 0 is the dominant simple eigenvalue of 
K with 1 = (1,1,..., 1) as a left eigenvector and a corresponding positive 
right (Perron) eigenvector, denoted by N, normalized so as to 1 • N = 1. 
The vector N = (N\, ■ ■ ■ ,N m ) is so-called the stable patch distribution; that 
is, the asymptotic, as t —> oo and disregarding demographic processes, dis¬ 
tribution of the population described by 183) among the groups. Under some 
technical assumptions 010, it can be proved that the total population 
n e = ni, e + ... + n mje converges as e 0 to the solution to 

OO 

d t n = — d a n — p*n, n(0,t) = J jd*(a)n(a,t)da, n(a, 0) = 1 • ri, (85) 

o 

where p* = p\Ni + .. .+^ m A m and /3* = f3\N\ + .. .+/3 m A r m are, respectively, 
the ‘aggregated’ mortality and fertility rates. We observe that, as required by 
the general paradigm, 1185)) does not display any explicit ‘space’ structure. 
However, a shadow of the original network persists in the model through the 
coefficients of the stable patch distribution vector N. More precisely, the coef¬ 
ficients of N give approximate fractions of the population residing on average 
in each patch and thus subject to patch specific mortality and birth processes. 
Therefore p* and ft* are aggregated death and birth coefficients which take 
into account that different fractions of the population die and give birth with 
different rates. 

At the same time we observe that despite being constructed as a 

microscopic extension of m, does not yield it back in the singular limit. 
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